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t—{ Abstract. We present a Ramsey- type atom interferometer operating with an 

optically trapped sample of 10 6 Bose-condensed 87 Rb atoms. The optical trap allows us 
to couple the \F = 1, mp — 0) ^ \F = 2, mp = 0) clock states using a single photon 
6.8 GHz microwave transition, while state selective readout is achieved with absorption 
imaging. Interference fringes with contrast approaching 100% are observed for short 
evolution times. We analyse the process of absorption imaging and show that it is 
possible to observe atom number variance directly, with a signal-to-noise ratio ten 
times better than the atomic projection noise limit on 10 6 condensate atoms. We 
discuss the technical and fundamental noise sources that limit our current system, and 
outline the improvements that can be made. Our results indicate that, with further 
experimental refinements, it will be possible to produce and measure the output of a 
sub-shot-noise limited, large atom number BEC-based interferometer. 
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1. Introduction 

Atom interferometers have evolved significantly over the last two decades, from systems 
with matter-based gratings and thermal beams, to systems based on ultra-cold atoms 
and optical or microwave transitions [1]. Atom interferometers are the basis of our 
time standard, in the form of laser-cooled Cs fountain clocks with accuracies of a few 
parts in 10 16 [2]. The sensitivity of inertial sensors based on atom interferometers is 
quickly approaching, and in some cases surpassing, state-of-the-art mechanical or optical 
systems [3]. At the same time a number of measurements have proven the competitive 
nature of atom interferometers in fundamental tests and measurements [4]. Perhaps 
more tantalisingly, it has been proposed that instruments based on atom interferometers 
can complement the LIGO gravitational wave observatory, allowing access to the 1 mHz 
to 10 Hz frequency band [5]. 

Laser-cooled atoms have become the source of choice for atom interferometry [1]. 
The very low centre of mass velocity and narrow energy width of laser-cooled samples 
allows for a more compact apparatus and better beam-splitting efficiency than hot 
thermal beam based systems. For optical interferometry, lasers are a very convenient 
source, offering high flux, low divergence and narrow linewidth. An analogous source 
for atoms is a Bose-Einstein condensate [6, 7, 8]. In addition to having a dramatically 
narrower velocity width than a laser cooled source, which facilitates high contrast, large 
momentum beam splitting, BECs offer the possibility of quadrature squeezing [9] as 
a path to improving the signal-to-noise in an interferometric measurement [10]. Very 
recently, an increase in interferometric sensitivity has been demonstrated via quadrature 
squeezing in small condensed samples [11, 12]. 

The possibilities offered by Bose-condensed sources have driven the development of 
magnetically trapped clocks and waveguide structures for interferometry with the aim 
of producing robust portable devices [13]. These devices typically operate with small 
atom numbers, N (on the order of 1000 — 10,000 atoms), for which projection noise 
fluctuations are a large fraction of the total atom number (y/~N /N ~ 1 % to 3 %) and 
thus are readily observable. Recently, our group has demonstrated a free space atom 
interferometer operating at the quantum projection noise limit with 10 4 atoms [14]. 

Ideally, one would like to use larger atom clouds for interferometric sensors, since the 
achievable phase sensitivity increases with number (A0 ~ l/y/N). Using condensates 
containing 10 6 atoms (readily created in many BEC machines) would immediately give 
more than an order of magnitude improvement in signal-to-noise over current trapped 
atom interferometers. The shot noise on such a large cloud is of order y/N /N ~ 0.1 %, 
and detecting at this level puts stringent requirements on all technical aspects of the 
experimental system. In particular, the interferometer beam-splitters must be extremely 
stable, and the state-selective measurement of atom number must be very precise. 
Understanding the technical limitations to these components via measurements and 
calculations is critical to progress on large condensate atom interferometry. 

This paper is divided into several sections: Section 2 introduces the sensitivity 
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of a Ramsey-type atom interferometer, including a theoretical analysis of quantum 
projection noise and typical experimental noise sources. Section 3 describes the 
apparatus, and presents results from our trapped atom interferometer. Section 4 
discusses the requirements for an absorption imaging system capable of measuring below 
the shot noise limit on a large atom cloud. We present results on shot noise limited single 
images as well as analysing the possibility of observing large atom number squeezing 
with absorption imaging. Section 5 proposes an interferometer with sensitivity enhanced 
beyond the standard quantum limit that could be implemented with the apparatus 
described in this paper. We conclude by predicting the achievable signal-to-noise ratio 
for clocks and inertial sensors based on condensed atom interferometry using the clock 
states. 




Figure 1. The Ramsey interferometer scheme. An atomic wavepacket is split into two 
components, allowed to evolve for a time T, and then recombined. The atoms can be 
coupled to a different internal state, remaining spatially overlapped, or can be coupled 
to another momentum state, so that the interferometer encloses an area, (a) A spatial 
Ramsey interferometer is sensitive to inertial effects, (b) The subject of this paper: a 
temporal Ramsey interferometer is sensitive to state dependent phase shifts. A n pulse 
allows reflection for the separated beam path interferometer and imposes a 'spin echo' 
effect for the trapped system, (c) The experimental setup. A single photon microwave 
transition drives internal state transitions in a BEC held in a crossed dipole trap. 



2. The sensitivity of an atom interferometer 

In a Ramsey atom interferometer, an input atomic wavepacket is split into two states 
which are allowed to evolve for a time T until being recombined (see Figure 1). In this 
section, we examine the sources of noise in a typical Ramsey interferometer, and present 
a rigorous analysis of how these affect the sensitivity. 
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2.1. Noise sources 

For a temporal Ramsey interferometer (Figure 1(b)), the noise may be categorised as 
follows: (i) quantum projection noise, due to the finite number of particles involved 
in the measurement; (ii) noise introduced by the beam-splitting process; (iii) noise 
affecting the phase evolution during the evolution time; and (iv) noise on the detection. 
We discuss each of these in turn. 

2.1.1. Quantum projection noise After the final beamsplitter of a n/2 — n/2 Ramsey 
sequence, the atoms exist in a superposition of two states |1) and |2). Upon 
measurement, each atom is projected onto either state |1) or state |2), with probabilities 
Pi and p 2 respectively. Due to the fundamentally random nature of this projection, 
the final number of atoms in each state will not necessarily equal the expected value, 
but will fluctuate around this value in any given measurement. This fluctuation is 
quantum projection noise. Experimentally, the quantity we analyze is the fractional 
population in state |2), p = p 2 = N 2 /N, where N = Ni + N 2 . In the absence 
of squeezing, the projection noise causes binomial fluctuations in this quantity of 
o~ p = a/p(1 — p) /N = | sin6\/(2y/N). Operating halfway up a fringe, where 9 — | 
and the change in population for a given phase shift 5p/5(j) is maximal, the population 
in each state is Ni — N 2 — N/2, and the projection noise is a p = \/(2y/~N). 

Quantum projection noise represents a fundamental limit to the sensitivity of an 
atom interferometer. The simplest way to reduce this noise is to increase the atom 
number, since a p oc 1/ y/N. This is the motivation for large atom number interferometry; 
an interferometer made with 10 6 atoms is intrinsically 30 times more sensitive than one 
using 10 3 atoms. As in optical interferometers, the projection noise may also be reduced 
beyond the limit imposed by quantum projection noise by quadrature squeezing [15]. 
This is discussed further in Section 5. 

2.1.2. Beam-splitter noise In a typical Ramsey interferometer, the first beamsplitter 
ideally produces an equal superposition \ip) = | 2) )/V% of the two states, known as 
a 7r/2-pulse. The second produces a superposition which depends on the phase difference 
(j> between the atomic transition and the coupling field accumulated during the evolution 
time T. As the phase of the second beamsplitter is scanned, the population in each 
state varies sinusoidally, producing Ramsey fringes. The precision of an interferometric 
measurement is limited by the repeatability of these superpositions - any fluctuations 
in the splitting translates into fluctuations in p at the output of the interferometer. The 
coupling between the two states can be realised via a single-photon process (for example, 
microwave coupling between the ground states of 87 Rb, as in this work) or a two-photon 
process (for example, Raman coupling using optical transitions to an excited state). 
In either case, noise in the power or duration of the beamsplitter pulses may limit the 
sensitivity of the interferometer. The frequency and intensity of the radiation coupling 
the states must therefore be kept stable, and the energy of the transition itself, which 



Optically trapped atom interferometry using the clock transition of large 87 Rb BECs 5 
may be affected by external electric or magnetic fields, should not vary 

2.1.3. Phase evolution noise By design, an interferometer is sensitive to differences 
in the phase accumulated by each state during the interaction time. To measure with 
high signal-to-noise, the effect being measured (acceleration, rotation, magnetic field, 
etc.) should exist on a small background of noise. However, any other process which 
affects the phase of each state differently will contribute to the noise on the output. 
Factors which affect the entire cloud uniformly, such as fluctuations in magnetic field, 
result in variations in the final population in each state and lead to an increase in noise. 
Inhomogeneous effects such as position-dependent shifts due to the trapping potential, 
local magnetic field inhomogeneity, and shifts due to the collective excitation of the 
spatial mode cause phase evolution during the time T to proceed at different rates 
across the trapped cloud. This leads to a classical dephasing and a decrease in fringe 
contrast. 

2.1.4. Detection noise In this and many other recent publications on atom 
interferometry, the number of atoms in each state is measured using absorption imaging. 
This technique utilises Beer's Law to extract the atomic column density from a sequence 
of bright field images. Because the information on atom density is extracted from bright 
images the classical and quantum noise of the light is coupled into the measurement. 
With careful experimental design, classical noise can be minimised, leaving only the 
shot noise of the light field contributing to the readout noise. We analyse this situation 
in detail in Section 4. 

2.1.5. Other contributions Even without performing interferometry, there are a 
number of experimental difficulties to consider in reliably holding and imaging a large 
two component atom cloud. We must maintain a non-zero magnetic bias field to 
suppress spin exchange collisions, and additionally shield the cloud from external RF 
and microwave fields that may redistribute the interferometric states. The switching of 
magnetic fields and absorption imaging must not cause measurement errors larger than 
one part in 10 3 . 

2.2. Sensitivity analysis 

In this section, a general expression is derived for the atomic state after a Ramsey 
interferometer sequence: identical pulses of arbitrary duration t which bound the free 
evolution time T. The noise in this quantity is then analysed by propagation of 
uncertainties. In the analysis we make extensive use of the Bloch sphere representation 
of a two-state quantum system (see Ref. [16]). 
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2.2.1. State evolution First we define a 7r/2 pulse as one which brings an initially 
vertical Bloch vector into the x-y plane, by setting the length of the pulse to 
cos' 1 (-e 2 ) 

tf = ^^' (1) 
where e = |A| /fi is the magnitude of the ratio between the detuning and the Rabi 
frequency and f2 R — v 7 ^ 2 + A 2 = Qy/l + e 2 is the total Rabi frequency f. 

After the general Ramsey interferometer sequence described above, the final 
longitudinal spin projection P z , i.e. the ^-component of the Bloch vector (related to 
the population in state |2) by P z = 1 — 2p), is given by 

P z = a 2 + (1 - a 2 ) sin (|A| T + f) , (2) 



where 



e 2 + cos (n R t) 

a = 



f(l + 2e 2 ) cos (fi R t) + 1\ Jl if sin(fi R t)<0 

n7r — tan , ; n~ ' 



e 2 + l 

2e v / TT^ 2 sin (fi R t) / ' | if sin (fi R t) > 

Using t = t n /2 from (1) in the general expression for the final longitudinal spin 
projection P z after the interferometer sequence, we have 

P z (t = t n/2 ) = sin [| A| T - sin' 1 (l - 2e 2 )] , (3) 

and so the Ramsey fringe first crosses P z = at the evolution time T = Tq given by 

|A| T = sin' 1 (1 -2e 2 ) . (4) 

The time T is that at which the interferometer is maximally sensitive to relative phase 
shifts. 



2.2.2. Effect of beamsplitter power fluctuations and resonance fluctuations Suppose 
the pulse sequence is calibrated for some detuning and Rabi frequency. We seek to find 
the change in the final longitudinal spin projection for fluctuations in the power 5P 
and detuning SA of the coupling field, while the duration of the pulses t = t n / 2 and 
evolution time T = T are kept fixed. As such, the change in the spin projection is 
found by evaluating § 

X In this analysis, < e < 1 since for e > 1 it is not possible to achieve a 50 — 50 beamsplitter. The pulse 
duration required for a ir/2 pulse is shortest on resonance, t* (A = 0) = 7r/(20), and approximately 
equal to this value when 6«1. 

§ To convert these derived uncertainties in the z-projection of the Bloch vector into uncertainties in 
the population of state |2), use 5P Z = 2Sp. 
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for arbitrary (but constant) t and T. Evaluating the expression (5) and then substituting 
t = £71-/2 and T = T gives 

-p) + ^ (e) U 



(6) 
(7) 



where 



9(e) = 



e [cos' 1 (-e 2 ) - VT^ 



:i+6 



2^/2 



sm 



l_2e 2 ) 2 [VT^ + e 2 cos' 1 (-e 2 )} 



+ 



For e < 1, this can be approximated as 



71 



- 1 



/7T\2 /<J A 



2^3/2 



(8) 
(9) 



+ LH 



(10) 



5P N 

+ \2eJ V ^ 

where the approximation is valid to within 5 % for e < 0.4. Each contribution to 5P Z is 
monotonic in e: fluctuations in P z due to power increase with increasing e and those due 
to resonance fluctuation decrease with increasing e. The minimum of 5P Z is achieved 
when d(5P z ) 2 /de = while holding SA/Q and 5P/P constant: 

2/n\2 / x / r, \ 2 



5' (e) 



where the approximation is valid to within 10 % for e < 0.32. The optimal detuning can 
thus be found: 



-opt 



5 A P 
~tt5P 



V 1 - 2/tt 



IA 



opt I 



1.664 



5A n 



5P 



which is valid to within 8 % for e < 0.32. 

The fluctuations in detuning SA depend primarily upon fluctuations in the applied 
magnetic field B and the applied pulse frequency The frequency difference between 
the |1,0) and the |2, 0) states in 87 Rb is found using the Breit-Rabi equation to be 
/ = /oVT+ B 2 x 2 where fo — 6.834 GHz is the hyperfine splitting at B — 0, and 
x = Hb(9j — 9i)j fo with \i& the Bohr magneton and gj and gi the Lande g-factors. 
Thus the fluctuation of the resonant frequency of the atoms will be 



5u TCS = 2n5f 

= 2nf Bx 2 
y/1 + B 2 X 2 

= k(B) 5B , 



5B 
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where k (B) is defined in the last line. The fluctuation in the detuning A = o; r es — cu ap p 
is then given by 

(SA) 2 = (5u TCS f + (5w ap p) 2 
= k 2 (5B) 2 + (5u app ) 2 . 

To minimise sensitivity to power fluctuations, a small detuning A during the ir/2- 
pulses is required. However, this necessitates a longer evolution time for the spin to 
precess by ~ ir/2, resulting in more sensitivity to resonance fluctuations. Using a 
small (ideally zero) detuning during the pulses and a different detuning Ai during the 
evolution time || results in minimal sensitivity to fluctuations. Repeating the above 
analysis for A = then gives 




where the second order power fluctuations have been included. In this way P z can be 
made insensitive to power fluctuations to first order. 

3. An optically trapped atom interferometer 

In this section we present our Ramsey trapped atom interferometer operating on the 
|F = 1, m F — 0) — > \F — 2, m F = 0) transition in Bose-condensed 87 Rb. 

3.1. Apparatus 

Our experimental apparatus for producing optically trapped BECs of 87 Rb and 85 Rb 
is described in detail elsewhere [17]. In brief, we first load a retro-reflected 3D 
magneto-optical trap (MOT) with 10 10 87 Rb atoms in 5 seconds from a 2D-MOT 
source. The atoms undergo polarization-gradient cooling for 20 ms before being loaded 
into a quadrupole magnetic trap (lOOG/cm field gradient, 1G = 10~ 4 T) in the 
\F = 1, m F = —1) state. The cloud is transported over 4 cm to a quadrupole- Ioffe 
configuration (QUIC) magnetic trap with a bias field of B = 3.4 G and trapping 
frequencies u z = 2n x 16 Hz and oj p = 2n x 156 Hz. Radiofrequency forced evaporation 
then reduces the temperature of the sample to 10 /iK over 15 s. A crossed-beam dipole 
trap from a 1090 nm, 20 W fibre laser is then abruptly switched on and a homogeneous 
bias magnetic field of 160 G is added along the z-axis. One beam of the dipole trap 
propagates along the weak axis of the QUIC trap with a waist of 150 /im, while the 
other is almost perpendicular and in the horizontal plane with a waist of 200 /im. 
Approximately 10 7 atoms are thus held in a hybrid trap in which the radial confinement 
is provided by the dipole trap and axial confinement by the magnetic field curvature of 
the QUIC trap. 

|| This can be achieved by changing the atomic resonance by e.g. applying a magnetic field, or by 
changing the frequency of the reference oscillator. 
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Here the atoms are cooled for a further 7 s by lowering the intensity of the dipole 
trap laser. During this time the QUIC magnetic field is also ramped to zero leaving the 
atoms in a pure crossed dipole trap (u} XjViZ = 2n x (50, 57, 28) Hz) with a homogeneous 
bias field of 160 G. Finally, the atoms are transferred to the mp = state via a 5 ms 
Landau-Zener rf sweep. The result is an optically confined Bose-Einstein condensate 
containing 2 x 10 6 |F = 1, m F = 0) atoms. For an interferometry run, the bias field can 
be set to any value between G and 200 G. In this work we ramp the bias field down 
to 4G. 

We create a trapped Ramsey interferometer by coupling the \F = 1, m F = 0} and 
\F = 2, mp = 0) hyperfine states with a single microwave photon at 6.835 GHz. The 
microwave signal is produced by a pulse-gated Rohde & Schwarz SMR20 microwave 
generator locked to a rubidium frequency standard (SRS FS725). The microwaves are 
delivered to the experiment via a simple quarter-wave dipole antenna, essentially a 1 cm 
piece of copper wire soldered to an SMA connector. The antenna is glued in place as 
close as possible to the atom cloud. With this setup, we are able to couple the clock 
states with maximum Rabi frequencies of around 5 kHz. 

A typical interferometry sequence is as follows. A 300 pis ir/2 pulse creates an equal 
superposition of the |F = 1, m F = 0) and \F = 2, mp = 0) states (hereafter referred to 
as the |1) and |2) states respectively). After an evolution time T, another tt/2 pulse is 
applied and the trap is switched off, allowing the atom cloud to expand freely. During 
this expansion, a strong magnetic field gradient is applied to spatially separate the states 
|1) and 1 2) using the second-order Zeeman shift. After 20 ms of expansion, a pulse of 
light resonant with the |F = 1) — >■ \F' — 2) transition is applied to optically pump the 
|1) atoms out of the lower ground state, before both states are imaged by absorption on 
the \F = 2, m F = 2) — > \F' = 3, m F i = 3) transition. 

Figure 2(a) shows Ramsey fringes in the fractional population of the |2) state after 
the second tt/2 pulse, recorded by scanning the detuning of the microwave pulses from 
the atomic resonance. Fringes are shown for several different values of the evolution 
time T. Measuring only the relative population pi = ^/(Aq + AT 2 ) allows run-to-run 
variations in total atom number to be normalized out. The total number of atoms in 
these measurements is Aq + N2 = (1.0 ± 0.1) x 10 6 . 

To assess the performance of our system relative to the shot noise limit, we operate 
the interferometer at mid-fringe with N = 5 x 10 5 atoms and an evolution time of 5 ms 
with no spin-echo pulse. We find that the run-to-run standard deviation is a factor 
of 15 higher than the projection noise limit for this atom number, corresponding to 
a fluctuation of 2.1 % in the splitting probability P2 and a measured interferometric 
sensitivity of 1.8mrad over 30 minutes. 

3.2. Fringe visibility 

For short evolution times T < 10 ms and total atom numbers ~ 10 6 , we observe high 
contrast interference fringes with a visibility approaching 100% (Figure 2(a)). As T 
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Figure 2. (a) Decreasing fringe visibility with increasing evolution time in a 10 6 
atom interferometer. The uncertainties are smaller than the data points. (6) Fringe 
visibility with (n/2 — n — tt/2) and without (it/2 — it/2) a 'spin echo' w pulse. Error 
bars represent primarily statistical uncertainties, systematic effects are dominated by 
number fluctuations which arc normalized out in the measurement of relative number. 



is increased to 30 ms, however, the visibility decays to below 20%. Some of this loss 
of coherence is due to deterministic effects such as inhomogeneous broadening of the 
transition frequency across the trapped cloud. It is possible to reverse this dephasing by 
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applying a 'spin-echo' n pulse between the two Ramsey pulses (Figure 1(b)). This has 
the effect of flipping the state vectors on the Bloch sphere, thus reversing their spread 
during the second half of the evolution time. Figure 2(b) shows the effect on fringe 
visibility of adding a spin-echo pulse to an interferometer with N — 1 x 10 6 atoms. The 
pulse clearly increases the coherence time, with the time taken for the fringe visibility 
to drop to 70 % increasing from 5 ms to 20 ms. 

We also contend with irreversible inelastic losses from the interferometer states. The 
magnetically sensitive states \F = 1, m F = ±1), \F = 2, m F = ±1} can be populated 
by spin-exhange collisions [18], as shown in Figure 3(a). The inelastic loss rates also 
differ between the two interferometer states. Figure 3(b) shows the populations in each 
state over 1 s in the optical trap. The decay of state |2) is markedly faster than of 
state |1) due to inelastic dipolar relaxation from the upper hyperfine level. These losses 
contribute to the loss in fringe visibility as the evolution time is increased, but are also 
compensated for by a spin echo pulse. This is because the assymetric losses result in 
the Bloch vector falling out of the x-y plane during the evolution time due to heavier 
loss from state |2). The echo acts to bring it above the x-y plane, where it then falls 
back into the plane prior to the second n/2 pulse. 

The remaining fringe decay is a combination of technical noise and irreversible 
dephasing due to atomic interactions [19]. One would expect the majority of technical 
noise (e.g. current noise in the magnetic field coil power supplies) to affect the entire 
cloud uniformly, which would cause an increase in the noise on the fringes but no decrease 
in visibility. To a good approximation, any remaining inhomogeneous technical noise 
should cause a decay in fringe visibility independent of atom number, while interaction- 
induced dephasing will affect large atom number condensates more due to their higher 
density (this will be discussed in detail below). For an interferometer with ^ = 5x 10 4 
atoms, we measure a decoherence time of over 800 ms, with the fringe visibility still at 
80 % after 200 ms. Recent work has also shown that a self-rephasing effect can lead to 
coherence times of many tens of seconds [20]. We have not yet observed this effect in 
our lab. 

Thus we conclude that, in our setup, interaction-induced dephasing is the limiting 
factor for achieving long coherence times with large atom numbers. This could be 
overcome by making the samples more dilute (using a weaker trap) or by reducing the 
interactions between atoms using a Feshbach resonance. Manipulating the scattering 
length has been demonstrated to give a large increase in interferometric contrast in 
Bloch oscillations [21, 22]. 

3.2.1. Effect of interactions To see how the interactions affect the fringe contrast in 
a simple 7r/2 — n/2 interferometer scheme, consider a two mode semiclassical model 
of the BEC system [23], where we ignore the dynamics of the spatial mode. In our 
interferometer the condensate is split into two spatially overlapping hyperfine states. 
Assuming there are no multimode excitations, the mean field wavefunction may be 
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Figure 3. Inelastic processes which lead to a loss in fringe contrast, (a) Integrated 
optical depth of a sample of N = (5.0 ± 0.2) x 10 5 atoms. The \F = 1, mp = ±1) and 
\F = 2, nip = ±1) states are populated by spin-exchange collisions, (b) Measured loss 
rates for the two interferometric states. The inset shows the interferometer states |1) 
and 1 2) after 0.5 s and 1 s hold times. The error bars represent the combined statistical 
and systematic uncertainties. 



written as 

*(r,*) = c 1 (^i(r)|l) + c 2 (*)V 2 (r)|2) , (12) 

where ^i(r) and ^(r) represent the normalised spatial mode functions of the the |1) 
and 1 2) components of the condensate. The number of atoms in each mode is given 
by | c x | 2 = N u |c 2 | 2 = N 2 , with iVi + N 2 = N, and / |\I/| 2 d 3 r = N, normalized to the 
total number of atoms. From the Gross-Pitaevski equation, we obtain the equations of 
motion for the coefficients C\ and c 2 : 
dc 

i ~dt = UlCl + ( 9llNl + 9^2)^ = u a c 1 (13a) 
= W 2C 2 + (922N2 + gi 2 N 1 )c 2 = to h c 2 , (136) 
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where hui^ are the linear contributions to the energy of each state, and 

9ij 'I'' I |^(r)|%(r)| 2 rf 3 r. (14) 



For large Bose-Einstein condensates, we can use the Thomas- Fermi approximation for 

where N is the total number of atoms, fi is the chemical potential, and Uu and V(r) 
denote the s-wave interaction constant and the external potential, respectively, for the 
|1) component. This gives 

f^2V5 fl5U n \ 2/5 fmcu 2 \ 3/5 



<* = %w[rr) (irj • < 16 > 

where u is the geometric mean of the trapping frequencies for each dimension. It should 
be noted that g^ is a function of the total number of atoms N, but not the number of 
atoms in each mode. 

The rate at which the relative phase between the two modes evolves is 
d<P 

— = u a -u b 
at 

= (wi - w 2 ) + (ffn - #12)^1 - (#22 - #12)^2 • (17) 

Up to this point the evolution of the relative phase is completely deterministic, and 
hence reversible. However, we now take into account that when the atoms are coupled by 
the beam splitter, there is an uncertainty in the number in each state. Assuming a 50/50 
beam splitter, we can represent this number uncertainty by setting Aq = N/2 + AN, 
and since the total number of atoms is a conserved quantity, N 2 = N/2 — AN. Inserting 
this into (17) gives 

^ = ("1 - <" 2 ) + ten - 012) (f + AN ^j ~ (922 - gi2) - Aiv) .(18) 

As fluctuations in number due to spin projection noise are binomial (Section 2.1), this 
uncertainty in the number difference will be AA^ ps \^N /2. Ignoring the deterministic 
part of the phase evolution, the uncertainty in the evolution rate is 

A(tu a - tu b ) = ^( 9ll - 2g 12 + g 22 )VN , (19) 

and so the phase uncertainty after time T is 

A0 = A(u a - u b )T = -(g u - 2g 12 + g 22 )TVN . (20) 

A similar analysis shows that fluctuations in the total number will lead to an additional 
phase diffusion of A<p = (gu — g 22 )TAN, where AA^ represents shot-to-shot fluctuations 
in the total number. 

Now, if we consider a spin echo interferometer [re/ 2 — n — ir/2), the phase 
diffusion due to fluctuations in the total number is removed, since the phase difference 
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accumulated before the spin echo pulse is cancelled by that accumulated after it. 
However, this does not apply to the fluctuations in the relative number difference, as 
can be shown by tracking the phase of C\ and c 2 . Before the n pulse, after evolving for 
a time t\ = At, the phase of each mode is 

M*i) = ("l + 9n + AN^j + g 12 - AN^j ^ At (21a) 

Mh)= (u 2 + g22(^-AN^+g 12 (^ + AN^At. (2lb) 

The spin echo pulse does two things: it reverses the populations such that the population 
of mode 2 becomes the population of mode 1, and it also performs the mapping c\ — > ic 2 , 
c 2 —> ic\. So, directly after the spin echo pulse, at time t 2 , we have 

Mb) = Mh) + \ (22a) 

<h{h) = Mh) + \ ■ (22&) 
Evolving by another time At, the phase of each mode becomes 

0i (* 3 ) = Mh) + (^i + gii - AN^j + g 12 + AN^j ^ At (23a) 

Mh) = <h(h) + ("2 + 922 + AN^j + g 12 - AN^j ^ At , (236) 

where we note that the populations have been reversed, i.e. the population in mode 1 
is now N/2 — AN. Taking the difference in the phases at t — t 3 gives 

A0 = 2 (t 3 )-0i(t 3 ) 

= 2(011 - 2#i 2 + # 22 )AiVAt 

= (<?n-2(?i 2 + £ 22 )v / ]VAt , (24) 

or, in the Thomas-Fermi limit, 

(an-2ai 2 + a 22 )2(15v/m/i 2 a; 3 aii) 2/5 At 

A0= a^ 7H TVV^' (25) 

where is the s-wave scattering length for collisions between state % and j. Thus, the 
phase diffusion arising from interactions and the uncertainty in the number difference 
cannot be removed with a spin echo, leading to an irreversible phase evolution in the 
system. The phase diffusion due to fluctuations in the total number can be removed 
by this technique. For the parameters used in this experiment, we calculate the phase 
diffusion rate to be ~ 50 mrad s _1 . However, it should be noted that this rate is small 
due to the favourable scattering properties of 87 Rb, since (a n — 2a i2 + a 22 ) / an ~ —0.02. 

The interactions in the condensate also cause an overall density dependent shift to 
the resonance frequency [24]. Here, we do not consider this effect to be a source of noise 
but rather a measurable systematic effect. 
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3.3. Miscibility 

The intra- and inter-species scattering lengths in the \F — 1, mj? = 0) and 
\F = 2, m F = 0) states are predicted by full coupled channel simulations [25] using 
the potentials of [26] to be a u = 100.9 ao, 012 = 98.9 ao, and 022 = 94.9 ao, where ao 
is the Bohr radius. The miscibility parameter /i = 011022/012 ~ 0.98 < 1 predicts that 
the states will be immiscible, since the cross-species repulsion exceeds the self-repulsion. 
Spatial dynamics driven by differences in scattering length are clearly observed in the 
magnetic states \F = 1, mp = —1) and \F = 2, mp = 1) after T ~ 10 ms, which causes 
a loss of interferometric contrast [27]. However, we do not observe any dynamics of 
the |F = 1, m F = 0) and \F — 2, m F = 0) states after up to Is evolution times (see 
inset of Figure 3(b)), despite being able to detect domain separation in other states. 
This suggests either that spatial dynamics are occurring (and contributing to the decay 
in fringe visibility as T is increased) but are corrupted by our imaging procedure and 
thus not observed, or that the scattering lengths in our system are different to the 
theoretically predicted values. If the clock states are indeed miscible, this would make 
them more favourable for interferometry since miscibility minimises mean field driven 
dephasing and allows good wavefunction overlap for the recombination pulse. 

3.4- Beam-splitter stability 

For small detuning A, the fluctuations in P z after a single n/2 beamsplitter are 
SP Z ps (2 — 7r/2)(A/fi) 5A/Q + (vr/4) SP/P. Operating on resonance thus makes the 
beamsplitter immune to all relative frequency fluctuations to first order and requires the 
minimum rf power for a given n/2 pulse time. This gives good suppression of frequency 
instabilities in the microwave source and fluctuations in the resonance frequency due to 
changes in the magnetic bias field. As a rough estimate, a relative transition frequency 
fluctuation between runs of 10 Hz, which could be caused by a 10 Hz uncertainty in 
the microwave frequency or a 2 111G uncertainty in the 4 G magnetic bias field, would 
cause only a 0.006% fluctuation in transition probability on resonance for a 50/50 
beam-splitter. This increases to a 0.1% fluctuation 100 Hz off- resonance. The effect 
of magnetic field noise is reduced further at lower fields, since the two mp = states 
are sensitive only to second order in the magnetic field. Our measured background 
fluctuations are 4 mG peak-to peak, predominantly at 50 Hz. Our home built bias field 
power supply has a measured rms current noise of 1 part in 10 6 . 

While frequency fluctuations may be suppressed by operating on resonance, the 
amplitude of the microwave pulse affects the transition probability on a single beam- 
splitter linearly. We have measured the power variations in our 300 /is microwave it/ 2 
pulses to be approximately 0.5%, which precludes the observation of atom shot noise 
from a single beamsplitter with more than 10 4 atoms. Indeed, we observe fluctuations 
a factor of 2.2 above atom shot noise in a n/2 pulse with 3 x 10 5 atoms (see Figure 4). 

When operated on resonance, however, the spin projection P z after a full 
interferometer sequence is insensitive to power fluctuations to first order (see Section 



Optically trapped atom interferometry using the clock transition of large 87 Rb BECs 16 




10 15 20 

Number of runs 



30 



Figure 4. Performance of a single beamsplitter, showing the relative population as a 
function of the number of runs. This data has fluctuations a factor of 2.2 above atom 
shot noise on the 3 x 10 5 atoms used. 



2.2.2). Thus the 0.5 % power variation translates to only a ~ 0.0015 % variation on the 
interferometer fringes, which is a factor of 30 below the projection noise limit on our 
10 6 atom condensates. 

4. Shot-noise limited absorption imaging 

Absorption imaging is one of the most commonly used methods of studying ultra-cold 
atom clouds, and is our primary tool for extracting information from our experiments. 
In order to directly observe squeezing and a consequent increase in the sensitivity of an 
interferometer, it is necessary that the noise introduced by the imaging system be well 
below the atomic shot noise fluctuations. It has been shown that using CCD cameras 
with high quantum efficiency (QE), it is possible to directly measure atom number 
fluctuations at and below the projection noise limit on small ensembles in a series of 
absorption images [28] . Here we investigate whether this is possible for large clouds. We 
show that using a CCD camera with a quantum efficiency of only 17 % it is possible to 
achieve a signal-to-noise ratio of nearly 10 4 , almost an order of magnitude above that 
required to observe atomic shot noise on a 10 6 atom condensate. This is enough to 
observe a 9 dB reduction in atomic shot noise. 

Absorption imaging of cold clouds is performed by analysing two consecutive CCD 
images, one containing a shadow of the cloud and the other a bright-field image in the 
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absence of the atom cloud. We denote the number of electrons in the array of wells on 
the camera by ei and e/. These two images are processed pixel-wise across the array to 
find the number of atoms N px imaged at each pixel by the following equation: 



where Co = 2-nP 2 /(3A 2 M 2 ) is the resonant absorption cross-section scaled by the pixel 
area P and magnification M, e sat is the electron count corresponding to the saturation 
intensity I sa t = 1.67mW/cm 2 , and L = (4 A 2 + r 2 )/T 2 accounts for a detuning 
A of the imaging light from the \F = 2, mp = 2) — > \F' — 3, = 3) transition, 
where T = 2n x 6.067 MHz is the natural linewidth. The relationship between the 
electron counts recorded at a pixel and the intensity I px incident on that pixel is 
e P rr = T\Prl px j '(foou), where r\ is the camera quantum efficiency, r is the exposure time, 
and hu is the energy per photon. The light intensity at the position of the atoms is 
further scaled by the magnification of the imaging system. 

Assuming that all classical noise can be removed from the image, either via careful 
experimental setup or by numerical post-processing, there are only two noise sources 
in an image of the interferometric states: projection noise on the relevant number and 
photon shot noise, which manifests as fluctuations in e px . Read noise, dark counts 
and digitisation error are all negligible in our setup. The atom number in each cloud 
fluctuates according to the projection noise a a = ^/Np (1 — p), assuming a coherent 
spin state. The electron number fluctuations are Poissonian if the laser light incident 
on the camera is Poissonian, and so a d = y/e^. This can be expressed as a fluctuation 
o^et in measured atom number of 



where the sum is taken over all pixels i in a region of interest. The total noise in 
the image is then a = a/o" 2 + <y\ et . To directly observe atom number fluctuations, the 
detector noise a^et must be well below the atom shot noise a a . 

The quantity o~det is completely determined by the parameters of the imaging 
system. It is therefore possible to predict in advance the best way to image a cloud 
of a certain size and number. We have performed numerical simulations of absorption 
images, by solving (26) for e/, to find the imaging parameters which minimise o~d e t- The 
simulation determines the number of electrons accumulated in each pixel of the CCD 
array during an absorption image of a typical condensate following ballistic expansion. 
Figure 5 plots the cumulative standard deviation as a function of the number of runs 
for an imaging system optimised to image a single Bose-Einstein condensate of iV = 10 6 
atoms released from our crossed dipole trap. With the atom number fluctuations a a 
set to zero, we find the lowest achievable noise in the images to be a = 113 atoms, 
implying a detection noise of o~det/N = 0.0113 %. The imaging parameters are: intensity 
I51 sa t on resonance, magnification 8, exposure time 100 /is, and expansion time 30 ms. 




(26) 




(27) 
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Notably, the simulation assumes a camera with a low QE of 17%, which can be a 
factor of twenty lower in cost compared to the high QE cameras used in previous work 
[28]. This level of detector noise would permit observation of atom shot noise on a 
10 6 atom condensate with a signal-to- noise ratio of 9, or an atomic squeezing-enhanced 
interferometric measurement with sensitivity 9dB higher than the projection limit. For 
comparison, with the same camera, an iV = 10 5 condensate can be imaged with a ratio of 
5 between atom and photon shot-noise, with a magnification of 2 at saturation intensity. 
This is the configuration used in the current work. 




10 20 30 40 50 60 70 80 90 100 

Run number 



Figure 5. Simulated detector shot noise, showing the cumulative standard deviation 
and in the inset the actual fluctuation in measured atom number from run to run in 
the simulation. 

In practice, there are several further considerations when attempting to observe 
atomic shot noise by absorption imaging. Firstly, it is imperative that all sources of 
classical noise be removed from the imaging system. In our experience, the required 
stability can be achieved by careful design (removing unnecessary optical elements, 
using high stability mounts, isolating the beam path from air currents, removing dust 
particles, etc.). Secondly, the absolute atom number must be accurately calibrated by 
an independent method. There are several ways to do this, including analysing the 
phase transition from thermal cloud to BEC, or measuring the atom shot noise as a 
function of number. However, the atom number need not be calibrated to the level of 
the shot noise fluctuations; typically the number will be known in absolute terms to the 
2 — 3 % level, giving negligible error on the predicted shot noise. Finally, the number of 
experimental runs performed must be sufficient for the measured standard deviation in 
the atom number to be an accurate representation of its asymptotic value. As can be 
seen in Figure 5, 20 — 30 runs is typically sufficient to have an excellent estimation of 
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the asymptotic variance. 

5. Interferometry beyond the standard quantum limit 

Recent theoretical schemes to produce squeezing in an atom laser [29], as well as 
the experimental demonstration of squeezing-enhanced sensitivity of trapped atom 
interferometry [11, 12], point the way to future advances in precision measurement. 
We will now consider the enhancement in interferometer sensitivity by using a simple 
on-axis twisting scheme. In this scheme, a double interferometer (ir/2 — n/2 — ir/2) is 
used. The first interferometer prepares the quantum squeezing, and the second uses the 
squeezed state to perform a measurement with a sensitivity greater than the projection 
noise limit. To investigate how this would enhance the sensitivity of our setup, we 
have used a simple theoretical model, which is described in [30]. The sensitivity in our 
estimate of the phase shift in the measurement interferometer is fundamentally limited 
by the projection noise in the number difference measurement. The uncertainty in the 
phase due to this effect is 

i d<t> i 

where as before N tota iP z = iV\ — N 2 at the interferometer output. For uncorrected 
atoms, A0 = l/y/N tota i [31]. 

Figure 6 shows the phase sensitivity as a function of the phase of the second 
interferometer. The quantum projection noise limit is shown for reference. The case 
shown gives greater than a factor of 2 enhancement in the phase sensitivity, which is 
equivalent to increasing the atom number by a factor of 4. As the effective interaction 
strength for 87 Rb is quite small, it takes a relatively long time (20 ms) to prepare this 
level of squeezing. The maximum squeezing achievable for this system occurs after ~ 50 
ms. However, we considered this evolution time unrealistic as effects such as particle loss 
and technical noise will start to become relevant over this time scale. If we employed a 
Feshbach resonance to enhance the effective interaction strength, it would be possible 
to achieve a phase sensitivity of A0 x \J~N = 0.16 with a 0.25 ms evolution time. Any 
further level of quantum squeezing would not give significantly further enhancement, as 
we would then become limited by the detection noise. 

6. Conclusions and Outlook 

We envision a sensor in which a squeezed state is created in trap and then the atom 
cloud is released into free fall to minimise dephasing. Consider, for example, an atom 
based gravimeter such as the laser cooled system operating in Paris [32]. A sample 
of atoms could be produced as in the current experiment and released. The falling 
cloud can be split, reflected and recombined with Raman beams. We calculate that for 
the large numbers of atoms used in our experiments, the bare, unmodified scattering 
lengths of the an, a 22 and a 12 states will lead to a significant and usable amount of 
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Figure 6. Theoretical calculation of normalized phase sensitivity as a function of the 
phase of the second interferometer. The solid line shows the case where an on-axis 
twisting scheme has been used to prepare quantum squeezing, the dotted line shows 
the standard quantum limit. The preparation time of the interferometer was chosen 
as 20 ms. All parameters are consistent with the experimental set up presented in this 
paper. 

squeezing in 10 — 20 ms. The squeezed state could then be used as the source for the 
atom gravimeter. In such a system we calculate a single shot accuracy of 10 _9 g for a 
free fall time of 100 ms, with a further factor to be gained from sub-shot noise statistics. 
Such a system, based on BEC, is currently under construction in our laboratory. 

In conclusion, we have presented results from a trapped, large atom number 
interferometer operating on the |F = 1, mp = 0} — > \F — 2, mp = 0) transition in 
87 Rb. We have demonstrated that simple absorption imaging with an inexpensive CCD 
camera can have a sensitivity of 10 dB beyond the standard quantum limit for 10 6 
atoms. We have discussed the principal sources of technical and fundamental noise in 
our system which must be overcome to achieve projection-noise-limited performance, 
and have suggested techniques for minimising these. 
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t-H Abstract. In our original paper [1], we presented results from a Ramsey atom 

interferometer operating with an optically trapped sample of up to 10 6 Bose-condensed 
87 Rb atoms in the rap = clock states. We were unable to observe projection noise 
fluctuations on the interferometer output, which we attribute to the stability of our 
microwave oscillator and background magnetic field. Numerical simulations of the 
Gross-Pitaevskii equations for our system show that dephasing due to spatial dynamics 
driven by interparticle interactions account for much of the observed decay in fringe 
visibility at long interrogation times. The simulations show good agreement with the 
experimental data when additional technical decoherence is accounted for, and suggest 
that the clock states are indeed immiscible. With smaller samples of 5 x 10 4 atoms, 
we observe a coherence time of r = 1-0^3 s - 
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1. Interferometer stability 

As noted in the original paper [1, §3.2.2], fluctuations in the frequency difference 
A between the microwave reference oscillator and the atomic resonance during the 
beamsplitter pulses affect the transition probability p after a full interferometer sequence 
to first order. In addition, fluctuations in A during the interrogation time T inevitably 
translate into an uncertainty in p, as this is precisely what the interferometer is designed 
to measure. The latter dominates when the beamsplitter pulses are short compared 
with the interrogation time, as in our system. Assuming ideal resonant tt/2 pulses, the 
transition probability after a complete Ramsey sequence is 

P=\[l + cos(TA)] , (1) 

where A here represents the detuning during the interrogation time. The uncertainty in 
p caused by detuning fluctuations SA is therefore 5p = \dp/dA\ 5A = ^\ sin(TA)| SA. 
At the midpoint of a Ramsey fringe, TA = (2n + l)7r/2 for integer n, so | sin(TA)| = 1 
and 5p = T5A/2. (In fact, this is an underestimate, as it neglects the effect of 
fluctuations during the pulses.) Comparing this to the projection noise limit a p = 
1/y/AN, we can place an upper bound on 8 A: 

5 A < --*-= . (2) 

Fluctuations in detuning larger than this will preclude projection-noise-limited 
performance in a Ramsey interferometer with N atoms and an interrogation time T. 
With T = 5 ms and N = 10 6 , the allowable fluctuation is 8 A < 0.03 Hz. This requires 
a relative frequency stability of 4 x 10~ 12 in the microwave source, and magnetic field 
variations below 7/zG at a bias field of 4G. 

Operating our interferometer at mid-fringe with N = 5 x 10 5 atoms and an 
interrogation time of T = 5 ms, we measured run-to-run fluctuations in the transition 
probability of a v = 0.011 after a complete n/2 — n/2 Ramsey sequence. This is a 
factor of 16 higher than the projection noise limit for this atom number. In light of 
the above discussion, we conclude that the performance is limited by the stability of 
the microwave source (specified at 2 x 10~ n over 1 second) and the magnetic bias field. 
Remarkably, the measured fluctuations in p are below what would be expected given 
our measured background field variations of 4mG, which should cause 5 A = 20 Hz and 
thus a substantially larger 5p = 0.3. This can perhaps be attributed to averaging of 
fast fluctuations over the interrogation time. Nonetheless, it is evident that achieving 
projection-noise-limited performance with large atom numbers will necessitate a high 
degree of stability in the magnetic environment, with low-noise current supplies and 
effective magnetic shielding. It will also demand a highly stable reference oscillator to 
effect coupling between the states. 
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2. Fringe visibility 

In the original paper, we speculated that the observed decay of fringe visibility with 
interrogation time was due to a combination of inelastic loss, interaction-induced 
dephasing and technical noise. Here we show that while spatial dynamics driven by the 
interparticle interactions can account for much of the visibility reduction, it is unlikely 
that inelastic losses have a significant effect. We also quantify the decoherence due to 
technical noise sources, both with and without spin-echo pulses. 



2.1. Differential loss 

In a simple two-mode model, the only agent that can cause a reduction in fringe visibility 
(aside from imperfect coupling pulses) is differential loss from the interferometer states. 
If one of the internal states experiences greater loss than the other, the superposition 
prepared by the first n/2 pulse - which initially has equal amplitudes in both states 
- will evolve to an unequal superposition by the end of the interrogation time. The 
interference at the final beamsplitter will then be incomplete, yielding fringes with 
decreased amplitude. 

Suppose that the losses are such that, at the end of the interrogation time, a 
fraction hi of the atoms in state \i) remain in the system. The character and time- 
dependence of these losses are immaterial; they may arise from two- or three-body 
processes, background gas collisions or scattering by photons from the dipole trapping 
laser. The state vector before the final beamsplitter is then 



\f2 \ -i\fk 2 e 

and the fraction of atoms found in state |2) at the output of the interferometer, 
normalised to the total number p\ + p 2 = \{k\ + k 2 ) is 

1 " 



P= 2 



1 + — COS . 

ki + k 2 



(4) 



giving a fringe visibility of 



V=2v^. 

h + k 2 w 

Thus, having twice as many atoms remaining in one state compared with the other 
reduces the fringe visibility by only 6% - even with 90% loss from one state (and none 
from the other) the fringe visibility would still be 57%. f 

In our system, the interferometer states |1) and |2) do suffer from differential 
loss, primarily due to inelastic spin-exchange collisions which populate the 

| This weak scaling of the fringe visibility is entirely analogous to the situation in optical heterodyning: 
here, the weaker 'signal' state is amplified by mixing with a stronger 'local oscillator'. Of course, 
while losses affect the fringe visibility only inasmuch as they differ between the two states, an overall 
decrease in atom number raises the quantum projection noise limit and can therefore also decrease the 
interferometer's sensitivity. 
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\F = 1, m F = ±1), \F = 2, mp = ±1) states as demonstrated in the original paper [1, 
Fig. 3]. However, the analysis above demonstrates that this should have a minimal 
effect on the fringe visibility - even after 1 second of evolution, when k\ ~ 0.45 and 
A; 2 ~ 0.18, the visibility should remain above 90%. Over a typical interrogation time 
of tens of milliseconds, the contribution of the measured differential loss to decreased 
fringe visibility should be negligible (< 0.1%): to explain the visibility V ~ 0.1 observed 
at T = 30 ms would require a loss rate 400 times larger in one state than the other. § 
We conclude that differential loss does not contribute significantly to fringe visibility 
reduction in our system. 

2.2. Spatial dynamics 

The external state of the atoms participating in the interferometer - and in particular 
its associated phase - can also affect the output. The external state of a two-component 
Bose-Einstein condensate may be written in spinor notation as 

•™ -($:!!)• (e) 

where the ipi(r,t) represent the spatial wavefunctions of the two interferometer states, 
which obey a pair of coupled Gross-Pitaevskii equations and are normalised such that 
l^ + I^T = N. These can be expressed in terms of the particle density rij(r, t) = 
and a phase 0j(f, t) asipi(r,t) = >Jni(r, t) e % ^ r,t \ The phase difference 4>2(r,t) — <j)i(r,t) 
between the spatial wavefunctions, which can vary in space as well as with time, cannot 
be distinguished from that accrued during the interrogation time due to the energy 
difference of the internal states. As noted in the original paper [1, §4.2], inhomogeneous 
broadening mechanisms such as the position-dependent ac-Stark shift due to the dipole 
trapping potential or local magnetic field inhomogeneities cause the phase evolution 
during the interrogation time to proceed at different rates across the trapped cloud, 
reducing the fringe visibility. Those that are constant in time can be reversed using a 
spin-echo pulse. 

Even without an externally-imposed phase gradient, however, inhomogeneous 
phase evolution (with accompanying spatial dynamics) can be elicited by interparticle 
interactions. To investigate the impact of spatial dynamics on fringe visibility, we 
simulate the evolution of the two-component system in the mean-field approximation. 
For this purpose, we assume ideal n/2 pulses that perform the mappings ipi — > 
(il>i—iip2)/V2, 4>2 — > (4>2—i4>i)/V2- The initial state is taken to be a pure |1) condensate 
occupying the ground state of the external potential; to make the computation tractable 
we assume cylindrical symmetry and trap frequencies u ZtP = 2n x {30, 55} Hz, such that 

§ It should be noted that this analysis neglects spatial dynamics stimulated by the loss of atoms from 
each state; particularly in the case of density-dependent two- and three-body loss, the removal of atoms 
changes the density distribution and can excite oscillations. However, this should represent only a small 
perturbation in the limit that the loss rate is small compared with the trap frequencies, as is the case 
^re (r^ < w/10). 
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oj matches that of our crossed dipole trap. After an initial n/2 pulse, the evolution 
of the two components is determined by numerically integrating the following coupled 
GPEs: 



ih 
ih 





- h 2 


dt 


2m 




1 h 2 


dt 


2m 



-V 2 + V + g 22 \ij 2 \ 2 + g 1 2\ip 1 \ 



4>i (7) 
^2 , (8) 



where the nonlinear terms proportional to with g^ = 4iih 2 aij/m describe the 
effect of interparticle interactions, characterised by the scattering lengths an = 100.9 ao, 
a 12 = 98.9 a , and a 22 = 94.9 a . After evolution of the states for a time T, another 
7r/2 pulse mapping is performed with phase (f) n / 2 relative to the first, and the number 
of atoms in each state determined from ipi and ip 2 . By scanning the phase (p w / 2 of the 
final beamsplitter, we can simulate Ramsey fringes and determine their visibility. 

Figure 1 shows the results of this simulation for an interferometer with iV = 10 6 
atoms and interrogation times up to T = 40 ms. For the scattering lengths given above, 
the miscibility parameter fi = a\\a 22 ja\ 2 ~ 0.98 < 1 predicts that the states will 
be immiscible, since the cross-species repulsion exceeds the self-repulsion. Indeed, it 
can be seen from the simulated density profiles that the two states spatially separate, 
developing a ball-and-shell structure over several tens of milliseconds. As they are 
initially overlapped and thus far from the two-component ground state, oscillations 
ensue. At times when the relative velocity of the two components is greatest, a large 
gradient manifests in the relative phase and the visibility of the Ramsey fringes is 
reduced. 

The simulated fringe visibility is plotted as a function of interrogation time in 
Figure 2 (solid line), together with the experimental data from our original paper 
[1, Fig. 2b]. A faster reduction in visibility (well-described by an exponential decay 
with time constant r = 14(1) ms) is observed in the experiment, which we ascribe to 
inhomogeneous broadening effects such as the ac-Stark shift of the trapping laser and 
technical noise. The simulated fringe visibility also exhibits a different shape to the 
exponential decay seen in the experimental data, and we do not observe the predicted 
revival in fringe visibility after T ~ 20 ms. Qualitative disagreement between the 
experiment and simulation could arise from the condensate not being in its ground state 
before the first pulse, which would inevitably affect the subsequent dynamics. Perhaps 
more importantly, the simulation is cylindrical symmetric; the asymmetry of our crossed 
dipole trap may tend to wash out the periodic rephasing of the spatial dynamics. || 

We can also use the simulation to examine the effect of a spin-echo pulse on 
the spatial dynamics. To do this, we assume an ideal it pulse which performs the 

|| Although the fringes are simulated by scanning the phase of the final beamsplitter rather than the 
detuning, we do not expect this (or the assumption of ideal and instantaneous ir/2 pulses) to significantly 
affect the outcome since the interrogation times are much longer than the duration of the beamsplitter 
pulses. This is verified by comparing the simulated fringe visibility with that obtained by solving the 
full GP equations including coupling terms and varying A - the results agree to within 2%. 
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Figure 1. The impact of spatial dynamics on fringe visibility in a Ramsey 
interferometer with N = 10 6 atoms. From top to bottom: two-dimensional density 
distributions, radial density profiles, relative phase profiles and Ramsey fringes 
obtained from a simulation of two-component condensate evolution using the Gross- 
Pitaevskii equation. The spatial structure oscillates around the two-component ground 
state - a ball-and-shell arrangement due to the large intraspecies repulsion. When the 
relative velocity of the two components is highest, the relative phase exhibits a large 
variation across the cloud, reducing the fringe visibility. 



transformation ipi — > —iip2, ip2 —iipi a t time t = T/2, and again simulate fringes by 
scanning the phase of the final beamsplitter. The results for a T = 20 ms interferometer 
sequence with a spin-echo pulse are shown in Figure 3; in this case the predicted 
fringe visibility is increased from 20% to 75% (c.f. Figure 1). The pulse clearly 
does not reverse the nonlinear spatial dynamics, but partially mitigates the phase 
gradient which develops during the first half of the interrogation time, resulting in 
higher visibility fringes. The cancellation is not perfect: the altered density distribution 
causes the phase evolution due to mean-field interaction to be different during the 
second half of the interrogation time, preventing complete rephasing. Nonetheless, the 
pulse does significantly increase the coherence time of the interferometer. The fringe 
visibility obtained from the simulation is overlaid with the experimental spin-echo data 
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Figure 2. Visibility of Ramsey fringes from the in-trap interferometer with (filled 
circles) and without (open circles) a 'spin-echo' n pulse. Note that the data for the 
spin-echo interferometers represent the inferred visibility, which gives a lower bound 
on the actual fringe visibility. The lines show the fringe visibility determined from GP 
simulations of the spatial dynamics as described in the text; solid - no spin-echo pulse, 
dotted - with spin-echo pulse. The results of the spin-echo simulation agrees with the 
experimental data when an additional decoherence e~ T l T is added (dashed line). 



in Figure 2 (dotted line). The simulated spin-echo visibility qualitatively resembles 
the experimental data much more strongly than without the n pulse, suggesting that 
reversible inhomogeneous broadening does indeed contribute significantly to the decay 
in visibility in that case. By assuming an additional decoherence not remedied by the 
spin-echo pulse, V(t) = Ve~ T l T with time constant r = (120 ± 20) ms, we find excellent 
quantitative agreement with our experimental data (dashed line). 

2.3. Small atom numbers 

One would expect irreversible inhomogeneous broadening (such as due to technical noise) 
to cause a decay in fringe visibility independent of atom number,^ while interaction- 
induced dephasing should affect large condensates more due to their higher density 
and correspondingly higher interaction energy. This hypothesis is substantiated by 
GP simulations which show that, although smaller two-component condensates still 
separate spatially, the timescale of the dynamics is much longer, resulting in lower 
relative velocities and thus reduced phase gradients. As a final test of the impact of 
interaction-induced spatial dynamics, therefore, we run the interferometer with a much 
smaller sample of N = 5 x 10 4 atoms. The resulting fringes are presented in Figure 4(a), 

Aside from a small effect due to the increased spatial extent of larger condensates; in the Thomas- 
Fermi limit, the condensate radius scales weakly with atom number as oc iV 2 / 5 . 
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Figure 3. Spatial dynamics during a Ramsey interferometer with a spin-echo n 
pulse. The pulse is applied after 10 ms of interrogation time, exchanging the density 
distributions of the two states (upper plots) and inverting their relative phase profile 
(lower plots). During the second half of the interrogation time, there is a partial 
rephasing of the wavefunctions. The dashed lines in the lower plots show the 
relative phase evolution in the absence of a spin-echo pulse, which by the end of the 
interrogation time spans most of 2tt across the cloud. 

displaying high visibility V ~ 90% for interrogation times up to T = 52.5 ms. We are 
unable to observe fringes at larger T due to crippling phase noise, which scatters the data 
randomly as seen at T = 102.5 ms. Nonetheless, this is not indicative of decoherence, 
as the output does not converge to p — 0.5 as with the larger condensates. For the 
purposes of understanding decoherence, we can take the range of the scattered data as 
an estimate of the fringe visibility. The visibility is plotted as a function of interrogation 
time in Figure 4(b). Again, the visibility decays faster than predicted by the simulation. 
An exponential fit of the form e~ T l T yields a coherence time of r = 0.30(5) s, a factor 
of 20 longer than for N = 10 6 atoms. 

Adding a spin-echo pulse to this interferometer allows us to estimate the true 
coherence limit of our system in a regime where interaction-induced dephasing is 
minimal. Figure 5 plots the fraction of atoms measured in state |1) at the output 
of the interferometer after a7r/2 — 7r — 7r/2 sequence as a function of interrogation time. 
Unexpectedly, this fraction drops below p = 0.5 after 100 ms - most of the atoms are 
observed in state |2) at T = 200 ms. This cannot be explained by a phase shift that is 
constant in time, since this would be cancelled by the spin-echo pulse. A likely candidate 
for this effect is a deterministic drift in the magnetic bias field. The difference in phase 
accumulated before and after the spin-echo pulse is 5(p = (A 2 — Ax)T/2, where Ai 
and A2 represent the average detuning during each half of the interrogation time. If A 
drifts at a constant rate v, then 8(f) = uT 2 /4: resulting in fringes whose spacing decreases 
with T. The solid line in Figure 5 is a fit of this form, including an exponential decay 
to account for decoherence. From the fit, we calculate a drift rate of approximately 
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Figure 4. (a) Ramsey fringes from an interferometer with N = 5 x 10 4 and varying 
interrogation time. Data points represent the measured fractional population in state 
1 2), p = N2/N, and the solid lines are sinusoidal fits. At long interrogation times, 
run-to-run resonance fluctuations randomise the phase thus increasing the noise on 
the fringes. Where sinusoidal fitting fails, we take the fringe visibility to be the range 
of the scattered data (dashed lines), (b) Fringe visibility as a function of interrogation 
time. The dotted line represents the visibility obtained from a GP simulation of the 
spatial dynamics. The solid line is a fit of the form e~ T / T , from which we extract a 
coherence time of r = (300 ± 50) ms. 



v = 2tc x 50 Hz/s, which at a bias field of 4G requires a drift in magnetic field of 
roughly 10 mG / s. This is consistent with the field drift measured in our system which 
we believe to be caused by warming of our bias field coils. Assuming that the measured 
points represent lower bounds on the fringe visibility at a given interrogation time, the 
reduction in visibility can be described by an exponential decay e~ T l T with a coherence 
time of r = 1.0±g;| s. 

2.4. Miscibility 

In our original paper, we reported that no spatial structure was observed in the two 
clouds at the output of the interferometer, and contemplated the possibility that the 
clock states were in fact miscible, in contrast to the prediction of the calculated scattering 
lengths. In light of the agreement between the experimental data and the simulations 
presented above, we believe it is likely that the states are indeed immiscible as predicted 
by coupled-channel calculations. The expansion of the condensate during time-of- 
flight, the Stern-Gerlach separation pulse and the repumping pulse applied before the 
measurement, along with the line-of-sight integration inherent in the absorption imaging 
itself, could all contribute to smoothing spatial inhomogeneities, concealing the structure 
developed during the interrogation time. As noted, states which are miscible are more 
favourable for interferometry, as only minimal evolution of the spatial distributions 
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Figure 5. Fraction of atoms measured in state |1) after a ir/2 — it — it/2 pulse sequence 
with N = 5 x 10 atoms. The population oscillates due to a systematic phase shift that 
is not cancelled by the spin-echo pulse, most likely due to a drift in the bias magnetic 
field changing the detuning during the interrogation time. An exponential fit to the 
envelope (dashed line) yields a coherence time of r = I.O+qJ s. 

occurs due to low-level breathing- mode oscillations excited by the first n/2 pulse. Figure 
6 shows the simulated fringe visibility for an iV = 10 6 atom interferometer with various 
values of the interspecies scattering length a± 2 . With an adjusted by only 10%, the fringe 
visibility does not drop below 90% over 100 ms of interrogation time. This suggests that 
an interspecies Feshbach resonance could be exploited to obviate interaction-induced 
dephasing and prevent the consequent decrease in interferometric sensitivity. 

[1] P.A. Altin, G. McDonald, D. Doring, J.E. Debs, T.H. Barter, J.D. Close, S.A. Hainc, T.M. Hanna, 
and R.P. Anderson. Optically trapped atom interferometry using the clock transition of large 
87 Rb Bose-Einstein condensates. New J. Phys., 13(6):065020, 2011. 
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Figure 6. The effect of interspecies interactions on fringe visibility in an in- 
trap Ramsey interferometer with N = 10 6 atoms, calculated from numerical GP 
simulations. The solid lines correspond to scattering lengths 012 where the two states 
are miscible; the dashed line shows the visibility for the background scattering length 
a 12 = 98.9 ao at which the states are immiscible. A reduction in a\i of only 10% is 
sufficient to maintain fringe contrast above 90% over 100 ms. 




